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Quantum chromodynamics (QCD) at finite temperature (T), baryon chemical potential (/ib) and 
isospin chemical potential (fii) is studied in the strong coupling limit on a lattice with staggered 
fermions. With the use of large dimensional expansion and the mean field approximation, we derive 
an effective action written in terms of the chiral condensate and pion condensate as a function of T, 
/xb and fix- The phase structure in the space of T and /xb is elucidated, and simple analytical formulas 
for the critical line of the chiral phase transition and the tricritical point are derived. The effects of 
a finite quark mass (m) and finite fii on the phase diagram are discussed. We also investigate the 
' phase structure in the space of T, fii and m, and clarify the correspondence between color SU(3) 

QCD with finite isospin density and color SU(2) QCD with finite baryon density. Comparisons of 
our results with those from recent Monte Carlo lattice simulations on finite density QCD are given. 
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CLi- I- INTRODUCTION 

<: 

■ To reveal the nature of matter under extreme conditions in quantum chromodynamics (QCD) is one of the most 
\ interesting and challenging problems in hadron physics. At present the Relativistic Heavy Ion Collider (RHIC) 
experiments at Brookhaven National Laboratory are running in order to create a new state of hot matter, the quark- 
gluon plasma [ij, |2( ■ On the other hand, various novel states of dense matter have been proposed at low temperature 
which are relevant to the physics of the interior of neutron star , such as the 3 P2 neutron superfluidity the pion £| 
and kaon [f| condensations and so on. Furthermore, if the baryon density reaches as much as ten times the nuclear 
density at the core of the neutron stars, quark matter in a color-superconducting state may be realized Q. 
1 One of the first principles methods to solve QCD is numerical simulation on the basis of the Monte Carlo method. 
y—{ \ This has been successfully applied to study the chiral and deconfinement phase transitions at finite temperature 
C*") ■ (T / 0) with zero baryon chemical potential (/ib = 0) 7}. On the other hand, lattice QCD simulations have 
an intrinsic difficulty at finite /xb due to the complex fermion determinant. Recently, new approaches have been 
proposed to attack the problem in SU C (3) QCD at finite /j,b @, such as the improved re weig hting method @, the 
£v' Taylor expansion method around /^b = [T(j and analytic continuation from imaginary //b jlll fl3 |. They have 
indeed given us some idea about the phase structure of QCD, e.g., the existence of the critical end point of the chiral 
phase transition and the determination of the slope of the critical line. They are, however, inevitably restricted either 
to simulations on a small lattice volume or to the low ^b region near the critical line. Alternative attempts have 
• • ■ also been made to study theories with positive fermion determinant, which include color SU C (2) QCD with finite /zb 
.£h ! 0113 and color SU C (3) QCD with finite isospin chemical potential pjlllpl. 

^ In this paper we consider the strong coupling limit of SU C (3) lattice QCD |l7t ITU lilj l2fl l2ll l22j | and make an 

extensive study of its phase structure with the use of the effective potential analytically calculated at finite T, /ib, /ii 
and the quark mass m. The phases obtained as functions of these variables are not easily accessible in Monte Carlo 
simulations and thus give us some insights into hot and dense QCD. Indeed, our previous study of strong coupling 
SU C (2) QCD with finite T, /iB and m show results which agree qualitatively well with those in numerical simulations 
[23|. Note that our approach has a closer connection to QCD than other complementary approaches based on the 
Nambu— Jona-Lasinio model |24|. the random matrix model [25| and the Schwinger-Dyson method p^.l2?|. 

This paper is organized as follows. In Sec.[nJ we consider color SU(7V C ) QCD with finite temperature and baryon 
density. Starting from the lattice action in the strong coupling limit with one species of staggered fermion (corre- 
sponding to Nf = 4 flavors in the continuum limit), we derive an effective free energy in terms of a scalar mode a and 
a pseudoscalar mode 7r for N c > 3 using a large dimensional expansion and the mean field approximation (Sec. Ill A"|l . 
The resultant free energy is simple enough that one can make analytical studies at least in the chiral limit m = with 
finite T and ^b- Secion ITTBl is devoted to such studies and useful formulas for the critical temperature and critical 
chemical potential for the chiral restoration are derived. Inpaticular, we can present an analytical expression for the 
temperature and baryon chemical potential of the tricritical point. In Sec. Ill CI we analyze the chiral condensate and 
baryon density as functions of /is and show the phase diagram in the plane of T and /zb for N c = 3, N{ = 4. We 
discuss the effect of finite quark mass m on the phase diagram of QCD and comparison with the results from recent 
lattice simulations are given here. 

In Sec. IIIII we consider the strong coupling lattice QCD with isospin chemical potential as well as T and /ib- In 
order to include we extend the formulation studied in Sec.[n]to two species of staggered fermion, which corresponds 
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to Nf = 8 flavor QCD in the continuum limit (Sec. IIII A"|) . Then we restrict ourselves to two particular cases: One is 
at finite T, finite /j,b and small /ij (Sec. 1111 B)) : the other is at finite T, finite \i\ and zero \ib (Sec. 1111 T^|) . For each 
case, we derive and analyze an analytical expression for the effective free energy and show the phase diagram in terms 
of T and /j,b or /ii for N c — 3, N{ — 8. In Sec. IIII Bl we discuss the effect of increasing Nf and that of finite fij on 
the phase diagram of QCD in the T-fiB plane. In Sec. IIII 01 comparison with the results from lattice simulations and 
discussion of a correspondence between QCD with finite /ii and SU C (2) QCD with finite /iB are given. In Appendixes 
^ and [51 we give some technical details for deriving the effective free energy. 

II. QCD (N c > 3, N t = 4) WITH FINITE BARYON DENSITY 

In this section, we consider strong coupling lattice QCD with finite baryon density. First of all, we review how to 
derive the effective free energy written in terms of a scalar mode a and a pseudoscalar mode 7r for N. > 3 for further 
extension in Sec. IIIII while the resultant free energy is almost the same as that studied in fl8L lisj . Then we 
analyze the chiral phase transition in the chiral limit at finite temperature and density, and derive analytical formulas 
for the second order critical line and the tricritical point. Finally we show the phase diagram in terms of T and fiB 
for N c = 3. 

A. Formulation of strong coupling lattice QCD 

In order to derive the effective free energy written in terms of the scalar and pseudoscalar modes, we start from 
the lattice action with one species of staggered fermion, which corresponds to the continuum QCD with four flavors. 
In the strong coupling limit, the gluonic part of the action vanishes because it is inversely proportional to the square 
of the gauge coupling constant g. Consequently, the lattice action in the strong coupling limit is given by only the 
fermionic part: 

S[U, x, x] = m J2 X(x)x(x) + 2 Voix) {x{x)e"U (x)x(x + 6) - %(x + 6)e-^U f (x) x {x)} 

X X 

1 . (1) 

+ 2 X^X^'^ {x(x)U j (x)x(x+j) - x(x + j)u](x)x(x)} . 

X j = l 

X stands for the quark field in the fundamental representation of the color SU(iV c ) group and is the SU(iV c ) 
valued gauge link variable, d represents the number of spatial directions and we use a notation x = (r, x) in which r 
(x) represents the temporal (spatial) coordinate. T]o(x) and T)j(x) inherent in the staggered formalism are defined as 

rjo(x) — 1, rjj(x) = (— l)^*=i x *- 1 , |U is the quark chemical potential, while the temperature is defined by T = (aNr)^ 1 
with a being the lattice spacing and N T being the number of temporal sites. We will write all the dimensionful 
quantities in units of a and will not write a explicitly. 

This lattice action has U(l)v x U(1)a symmetry in the chiral limit m = 0, which is a remnant of the four flavor 
chiral symmetry in the continuum theory, defined by U(l)v : x( x ) l— * e x( x ), x( x ) l— * x( x ) e ~ Wv an( i U(1)a : x( x ) l— * 
e le ( x ' eA x(x), x(x) l— ¥ x(x)e 1E ( x ' 6a . Here e{x) is given by e(x) — (— l^^o 2 ^ which plays a similar role to 75 in the 
continuum theory. U(1)a will be explicitly broken by the introduction of a finite quark mass or spontaneously broken 
by condensation of (xx)- Note that the staggered fermion's U(1)a symmetry should not be confused with the axial 
U(l) symmetry in the continuum theory, which is broken by the quantum effect. Now by using the lattice action 
Eq. m, we can write the partition function as 

Z = J V^xMUomU^e- 3 . (2) 

In the succeeding part of this subsection, we perform the path integrals by use of the large dimensional (1/d) expansion 
and the mean field approximation to derive the effective free energy. 

First, we perform the integration over spatial gauge link variable Uj using the formulas for the SU(iV c ) group 
integration. Keeping only the lowest order term of the Taylor expansion of e~ s , which corresponds to the leading 
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order term of the 1/d expansion [23, we can calculate the group integral as follows: 



T>[Uj} exp 



1 d 

2 ^2^2vj(x) [x(x)Uj(x)x(x+j) - x[x + j)U){x)x{x)} 



x j=l 



exp 



Y,M(x)V M (x,y)M(y) 



(3) 



where M{x) and Vm(x, y) are the mesonic composite and its propagator in the spatial dimensions, defined, respectively, 
by 



M{x) = ±-5 ab x a {x)x\x), 



V M (x,y) = ^Y,{ 6 : 



^ + 5 

y,x+j y,x-j 



(4) 



The resultant term in the action describes the nearest neighbor interaction of the mesonic composite M{x). Residual 
terms are in higher order of the 1/d expansion and vanish in the large dimensional limit. Note that Eq. © is correct 
as long as iV c > 3. For N c = 2 the baryonic composite (diquark field) gives the same contribution as the mesonic one 
in the 1/d expansion. The finite density SU C (2) QCD with a diquark field was studied in a previous paper p3 |. 

Next, we linearize the four-fermion term in Eq. @ by introduction of the auxiliary field a(x). This is accomplished 
by the standard Gaussian technique; 



exp 



J2M(x)V M (x,y)M(y) 



J V[a] exp 



{v(x)V M (x, y)a(y) + 2a{x)V M (x, y)M(y)} 



x,y 



(5) 



From the above transformation, we can show the relation between the vacuum expectation values of the auxiliary 
field and the mesonic composite: (c(x)) = — (M(x)). In order to keep the action invariant under the U(l)v and 
U(1)a transformations, we should assign the transformation for the auxiliary field as U(l)v : &(x) <— > cr(x) and 
U(1) A : <r(x) ^ e 2i£( - x}0A <r{x). 

In the staggered fermion formalism, one species of staggered fermion at various sites is responsible for the degree 
of freedom of flavors and spins. Therefore the auxiliary field cr(x) is responsible for both the scalar mode a and the 
pseudoscalar mode ir psj . In order to derive the effective free energy written using the scalar and pseudoscalar modes, 
we introduce a and it as <r(x) = a + ie(x)TT |29j. Then we can find a correspondence between these new values a, tt 
and the vacuum expectation values of the mesonic composites: a = — (x° '{x)x a {x)) /N c , n = (i \ a (x)e(x)x a (x)) /N c . 
These equations tell us that a and n are even and odd under lattice parity, respectively [l3l l30l | , and therefore a and 
7r correspond to the spatially uniform condensates of the scalar mode (chiral condensate) and the pseudoscalar mode 
(pion condensate). They transform under U(1)a transformation as 



/cos2#a — sin2#A\ (o 



cos 26 a J yvr 



U(1) A : 

The partition function is then written as 

Z = J V[U Q ]V[x,x] e- 5 '^-****] , 



(6) 



(7) 



with 



?'[Z7 , x, XI °M = rri X {x)x{x) + i [x{x)^U {x) X {x + 6) - X (x + 6) C -^U f (x) X {x)} 



N c d , , o \ N c d . , . , 
+ ^L. (cr +n 2 ) + ^-{a- ie{x)n} M(x) 



(8) 



Note that this action contains hopping terms of Xi X only in the temporal direction after the 1/d expansion and the 
mean field approximation. 

In order to complete the remaining integrals, we adopt the Polyakov gauge |l7j| in which Uq(t, x) is diagonal and 
independent of r: 



N,. 



U (t,x) = diag 



j(t>i(x)/NT 



,e 



i0iY c (5)/AT T 



with 4>a{x) = . 



(9) 
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Also we make a partial Fourier transformation for the quark fields; 



x(t, x) 



fWr 



X(r, x) 



fWr 



E 1 



-ik n r 



X(n,x) 



2tt 



(10) 



Substituting Eqs. J§J) and (IX Ul) into the action @ and taking the summation over r, the Grassmann integration over 
the the quark fields x an d X results in the following determinant: 



N c N T /2 

nnn 

x a—1 n—1 



sin kn 



M 



7-2 



(11) 



with k n = k n + 4> a {x)/N T — i/i. M denotes the dynamical quark mass defined by M = m + (d/2)a. 

The product with respect to the Matsubara frequencies n can be performed using a technique similar to that in 
the calculation of the free energy in finite temperature continuum field theory. The details of the calculation are 
given in Appendix 1X1 The result turns out to have a quite simple form: 2cosh[N T E] + 2cos[(j> a (x) — iN r /j] with 
one-dimensional quark excitation energy E[a, tt] = arcsinh[-\/A/ 2 + (d/2) 2 ir 2 ]. 

Finally we can perform the integration with respect to Uq by use of the formula in Appendix [B] for the Polyakov 
gauge. Then the integration gives, up to an irrelevant factor, 



/ n ri < 2 c ° sh + 2 co£ #a(^) - iN rtA} = n fe ^t p n+i _j 



(12) 



where Pq = 2 cosh [N T E], P±\ = cosh [N T fj] ± sinh [N T fj] = c ±Nrt * and i^„|>3 = 0. The determinant is to be taken 
with respect to i, j — 1,2,..., N c . The determinants have nonvanishing values only for n = 0, ±1 and their summation 
can be calculated exactly for general JV C . First, for n — 0, det Pi-j is expressed in the form of an N c x N c matrix as 

•• 



detPi_j = 



P P-l 

P+i Po '■■ 

'•■ '■■ P-! 

: P+i Po 



(13) 



Making a recursion formula for it and using the fact that P_iP+i = 1, we can obtain the simple solution det Pi— j 
sinh [(N c + l)N T E] / sinh [N T E] [lj|. Then the calculation for n = ±1 is rather easy and results in detP~i+i-j 

2 cosh [N c N T fi\. As a result, the effective free energy is given as follows; 



det P 1+i -j = P 



P 



N c d 



P cff [a, tt; T, Mb ] = (- log Z) / (EJ = ^ ^ 



-Tlog<^2cosh[/a B /T] 



sinh[(A^ c + l)£:/r]\ 



sinh [E/T] 



(14) 



Here we have defined the baryon chemical potential /jb as hb = N c fi. Although the temperature T = 1/N T takes 
discrete values, the effective free energy F e s[l/N T ] can be uniquely extended to a function of continuous T in terms 
of an analytic continuation since F e g[l/N T ] is defined on an infinite sequence of points. Hereafter we consider T as a 
continuous variable. This effective free energy in Eq. I|14f> corresponds to the Ni = 4 flavor QCD in the continuum 
limit. 



B. Analytical properties in the chiral limit 



The effective free energy obtained in the previous subsection is simple enough to make analytical studies in the 
chiral limit. This is one of the main advantages of our approach. Such analytical studies on the chiral phase transition 
are useful in understanding the phase structure of strong coupling lattice QCD, which will be presented in the next 
subsection. In the chiral limit m = 0, the quark excitation energy reduces to E[a : n] — arcsinh [(d/2) \J a 2 + ir 2 \ . 
Therefore the effective free energy in Eq. i|14|) is a function of only ^/a 2 + tt 2 . Of course this symmetry comes 
from the chiral symmetry U(1)a of the original lattice action. Since we can arbitrarily choose the direction of the 
condensate, we take a ^ and tt = 0. 
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Chiral restoration at finite temperature 

As we will confirm numerically in the next subsection, the effective free energy exhibits a second order phase 
transition at finite temperature. Therefore we can expand it in terms of the order parameter a near the critical point. 
Expansion of the effective free energy up to second order of the chiral condensate a gives 

F cS [a,0] = — a - 2mNc + 1 + 2cQsh[fiB/T] f +0{a). (15) 

As long as the coefficient of a 4 is positive, the condition that the coefficient of a 2 is zero gives the critical chemical 
potential for the second order phase transition as a function of temperature: 



fi^(T) = Tarccosh 



(JV C + l){d(N c + 2) - 6T} 
12T 



(16) 



The critical temperature at zero chemical potential is given by solving the equation /ig'(T) = and turns out to be 
r c (0) = d(JV c + l)(JV c + 2)/[6(JV c + 3)]. 

When the coefficient of a 4 becomes zero, the phase transition becomes of first order. An analytical expression for 
the tricritical point (T tr i, /Lttri) can be obtained as a solution of the coupled equations of Eq. (|TB|> and d 4 F e g/da 4 = 0, 
which results in 

rp _ V225iV c 2 + 20rf 2 (3A 2 +6N c -4)- 15N C _ 

J tri - 2Qrf ' Mtri — Mb Utrij • (i-i) 

The effective free energy exhibits a second order chiral phase transition for T > Ttri and it becomes of first order 
for T < Ttri- The existence of the tricritical point is consistent with the results of other anal ytic al approaches using 
the Nambu— Jona-Lasinio model |24|, random matrix model [25j. Schwinger-Dyson equation |26j and other methods 



Chiral restoration at finite density 

An analytical study of the first order chiral phase transition at finite temperature for T < T tr i is rather involved. 
However, the effective free energy reduces to a simpler form at T = 0: 



F e g[T = 0] = ^j-o- 2 - max {/x B , N c E[a]} , E[a] = arcsinh 



(18) 



It is easy to study the first order phase transition analytically in this case. 

The effective free energy has two local minima as a function of a: One is F e g = — at er = 0; the other is 
F c g = N c d<jQ /4 — N c E[ao] < at a = o-q, where ctq is the solution of the chiral gap equation dF e f[/d(T = 0, 



= 1 . 



therefore ctq = 



2vT + d 2 - 2 
J 2 



(19) 



As hb becomes larger, the global minimum changes from A c dCT 2 /4 — N c E[uq\ to —fiB at some value of the chemical 
potential. This is the critical chemical potential, given by ^™(T = 0) = N c E[ao] — N c da 2 /4. At this critical point, 
the order parameter a changes discontinuously as 



do when < ^b' 
when ^b > Mb" 



(20) 



This is a first order phase transition. We can also easily calculate the baryon density pb at T = 0, and it has a 
discontinuity associated with that of the chiral condensate: 



PB 



dF cS when ^ B < fjg 1 
df^B \ 1 when fi B > Mb' 



(21) 
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r=i.o r=o.8 




Baryon Chemical Potential jU B Baryon Chemical Potential jU B 

FIG. 1: Chiral condensate a and baryon density pa as functions of the baryon chemical potential pb for iV c = 3, iVf = 4. In 
the left panel, they are plotted just above the tricritical temperature (Tm = 0.866) at T — 1.0 for m — (the solid line) and 
m = 0.1 (the dashed one). In the right panel, they are plotted just below T tr i at T = 0.8 for m — (the solid line) and m = 0.1 
(the dashed one). 



At the critical chemical potential, pb changes from the empty density to the saturated density 1 |2fJ . Note that one 
baryon (= N c quarks) at one lattice site, that is pb = 1, is the maximally allowed configuration by the Fermi statistics 
for one species of staggered fermion. Saturation at large chemical potential on the lattice has also been observed in 
lattice simulations and strong coupling analyses of SU C (2) QCD |14ll23|. 

This chiral phase transition with finite pb at T = is not associated with the deconfinement phase transition to 
quark matter. Equations (|20|) and (|21[1 show that the chiral restoration coincides with the transition from vacuum 
to saturated nuclear matter, although they are in general separated p5| . Including the contribution of dynamical 
baryons, which is suppressed in the leading order of the 1/d expansion, will introduce the Fermi surface of baryons 
and hence separate the two phase transitions. 

C. The phase structure with finite quark mass 

In this subsection, setting the number of colors and spatial dimensions to 3, we examine numerically the nature of 
the chiral phase transition and the phase diagram of strong coupling lattice QCD for iV c = 3, Nf = 4. 

Chiral condensate and baryon density 

By minimizing the effective free energy F e g for N c = 3, d = 3 in terms of the order parameter a, we determine the 
chiral condensate as a function of T and pb- The baryon density pb is also calculated by —dF c ff/dpB- In Fig. ^ the 
results with zero and finite m are shown as functions of pb- 

In the left panel of Fig.Q] a and pb are shown for T = 1.0 just above T trl — 0.866. In the chiral limit they show a 
second order phase transition at the critical chemical potential given by Eq. 116|) , while introduction of an infinitesimal 
m makes the transition a smooth crossover. 

cr and pb are shown for T = 0.8 just below T t ri in the right panel of Fig.^ As discussed in the previous subsection, 
they show jumps at the critical chemical potential in the chiral limit. This first order phase transition is weakened by 
the introduction of m and finally becomes a crossover for large quark masses. 

The phase diagram 

Now we show in Fig. [5] the phase diagram of the strong coupling lattice QCD with N c — 3, Nf — 4 in the plane of 
T and pb- In the chiral limit, the chiral condensate a has nonvanishing value at low T and small pb, which means 
a spontaneous breakdown of the chiral symmetry Ua(1) by the condensation a ^ 0. The symmetry broken phase 
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FIG. 2: Phase diagram of strong coupling lattice QCD with iV = 3, Nf = 4 in the T-/^b plane. The solid (dashed) line denotes 
the phase transition line for the first (second) order phase transition in the chiral limit. The black circle represents the tricritical 
point (TCP) where the nature of the phase transition changes from first order to second. The first order line and its critical 
end point (CEP) for m — 0.4 are also shown. The dotted line denotes the flow of the critical end points when the quark mass is 
enlarged from 0. White circles are put on the critical end points for m = 0.1 (the left one), m = 0.4 (the middle) and m = 0.8 
(the right one) in the left panel. The vicinity of the TCP is enlarged and the first order line and CEP for m = 0.01 are shown 
in the right panel. 



is separated by a phase transition from the chiral symmetric phase; the phase transition is of second order in the 
high temperature region while it is of first order in the low temperature region. The second order critical line merges 
smoothly into the first order one at the tricritical point (T tr i,/Xtri) = (0.866, 1.73). 

In contradiction to our result, the chiral phase transition of four-flavor QCD in the chiral limit is thought to be of 
first order 0, • This difference comes from the fact that the symmetry of the staggered fermion action is limited 
to U(l) x U(l) in the strong coupling limit and we have employed the mean field approximation. 1 However, it is 
sensible to consider that the second order critical line in our result corresponds to the phase transition line of chiral 
restoration in QCD. 

As we can see from Fig. [21 the critical temperature falls much more ra pidly as a function of /j,b in contrast to the 
results from recent Monte Carlo simulations at small chemical potential [tj, U£l GjjLL2| • In order to make the discussion 
quantitative, we calculate the slope of the critical temperature near zero chemical potential. From Eq. I|16f> . 

= i _ 1.5-^ - 3.375-^ + O ( (u/T c ) 6 ) (22) 
T c (0) T c (0) 2 T c (0) 4 V vw c> J v ' 

with T c (0) = 5/3. Here we have written Eq. (12211 in terms of the quark chemical potential /i = ^b/N c in order to make 
the comparison to the lattice data easy. The value of the coefficient of the second term is two orders of magnitude 
larger than that calculated in the lattice simulation with N{ = 4 for m = 0.05 |12j . This large difference in the slope 
of the chiral phase transition may be understood due to the effect of m = and that of g — > oo. Our expression 122|) 
is derived in the chiral limit, while all of the lattice simulations are performed with finite quark mass. As reported 
in |34|. the chiral transition line has significant quark mass dependence in the lattice simulations: The slope of T c (/j,) 
suffers a sharp increase on decreasing m, from 0.025 (m = 0.1) to 0.114 (m = 0.005) in Nf = 3 QCD. In addition to 
that, as we can find in [l9|. the introduction of the next to leading order term of l/g 2 and 1/d in strong coupling 
lattice QCD decreases the critical temperature, and the amount of decrease at zero chemical potential is more than 
that at finite chemical potential. Therefore the slope of the critical temperature will be made softer by the corrections 
for the strong coupling and large dimensional expansions. Combining the above two effects, the difference of the slope 
between our result and lattice simulations could become reasonably small; this needs further study. 

Next we consider the effect of quark mass on the phase diagram. Introduction of finite m washes out the second 
order phase transition as shown in the left panel of Fig.^ As a result, the first order phase transition line terminates 



1 As one gets closer to the continuum limit, the symmetry will turn into the usual four-flavor chiral symmetry. Then the chiral fluctuations 
may induce a first order phase transition, as they do in the corresponding linear sigma model l33l . 
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TABLE I: Tricritical point (Ttri,ptri) for m = and critical end points (T on d, p-cnd) for various small quark masses. The 
deviations of the critical end points from the tricritical point are also shown. 



m 


rp 

J- end 


/^end 


1,-tend — -ttri^/Jtri 


(Mend — Ptri)/ptri 





0.866 


1.73 






0.001 


0.823 


1.73 


-5.03 % 


-0.541 % 


0.01 


0.764 


1.75 


-11.8% 


0.599 % 


0.05 


0.690 


1.85 


-20.4 % 


6.37 % 


0.1 


0.646 


1.96 


-25.5 % 


12.9% 


0.2 


0.590 


2.16 


-31.9% 


24.7% 



at a second order critical end point (T en d, Mend). As long as the quark mass is very small, m < 0.001, the critical end 
point flows almost along the tangent line at the tricritical point as in the right panel of Fig. [21 Then for the larger 
quark mass m > 0.001, Tend decreases and p cn d increases. Also, the first order phase transition line shifts in the 
direction of large pb associated with this change. 

Strong coupling lattice QCD may give us some indication for lattice simulations about the effect of m on the phase 
diagram of QCD. As shown in Tabled a finite quark mass m — 0.1 decreases T en d by 25% and increases Mend by 12% 
from the tricritical point (ltd, Mtri)- We can find in the same way that even a quark mass as small as m = 0.001 can 
move the critical end point by 5%. This indicates the relatively large dependence of the phase diagram of QCD on 
the quark mass. 

Our phase structure is consistent with those from other analytical approaches 0> HE HE E3- EUL E3 > except for 
one main difference j2^|. This is the fact that the gradient of the phase transition line near the tricritical point or 
critical end point is positive. Consequently, the critical end point flows in the direction of smaller T and smaller pb 
as a function of small quark mass. The positive gradient of the chiral phase transition line near the tricritical point 
can be understood as follows. The gradient dT/dp is related to discontinuities in baryon density Ap and in entropy 
density As across the phase boundary through the generalized Clapeyron-Clausius relation plf 

dT Ap 

dp As y J 

As we can see in the right panel of Fig. ^ the baryon density in the chiral restored phase p+ is larger than that 
in the chiral broken phase p_ just below the tricritical point (T = 0.8, pe = 1-7). We can also calculate the 
entropy density from the effective free energy (|14|> showing the result that As < 0. Ap = p + — p_ > holds 
because baryons are massless in the chiral restored phase and we can fill more baryons under fixed pe- Since 
we are in the leading order of the 1/d expansion, all baryons are static. Then we can estimate the entropy as 
s± ~ — p±logp± — (1 — p±)log(l — p±). Because 1 > p + > p_ > 0.5, the discontinuities in the entropy density is 
negative As = s + — s_ < 0. Therefore the gradient of the chiral phase transition line just below the tricritical point 
turns out to be positive, dT/dp = — Ap/As > 0. 



III. QCD (No > 3, Nf = 8) WITH FINITE ISOSPIN DENSITY 



In this section, we consider strong coupling lattice QCD with finite isospin density. First of all, we derive the 
analytical expression for the effective free energy in two particular cases. One is at finite T, finite pe and small pi 
with chiral condensates. The other is at finite T, finite pi and zero pe with chiral and pion condensates. Then we 
analyze the effective free energy and show the phase diagrams in terms of T and pe or pi for N c ~ 3, Nt — 8. 



A. Strong coupling lattice QCD with isospin chemical potential 



In order to investigate the effect of the isospin chemical potential pi on the phase structure of QCD, we have 
to extend the lattice action studied in the previous section to that with two species of staggered fermion having 
degenerate masses. The lattice action with two species of staggered fermion corresponds to the eight-flavor QCD in 
the continuum limit. Then for m = and pi = 0, it has U(2)y x U(2)a symmetry, as a remnant of eight-flavor chiral 
symmetry, defined by x(x) >-> e ilev+e{x)0A] ' T x(x) and x(x) ^ x(x)e~ ilev ^ £{x)eA] " r with r £ U(2). 
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After integrating out the spatial link variable Uj in the leading order of 1 / d expansion and introducing the auxiliary 
fields, we have the following partition function: 



V{U ]V[ Xl x}V{<j} e 



-S'[U , x ,x;<t] 



(24) 



with 



a—u,d x.y L 

+ 2l ^2[<yp a (x)V M (x,y)(7 a p(y) + 2<j,3 a (x)VM(x,y)M a f3(y)] . 

a,/3=u,d x.y 



x a (y) 



(25) 



Here the subscripts a, (3 represent the species of staggered fermion, taking "u" (for up quarks) and "d" (for down 
quarks). /i u , fid are the quark chemical potentials for each species of staggered fermion. o a p(x) are auxiliary fields for 
the mesonic composites M a p{x) = Xa( x )xp(x)/N c , and these vacuum expectation values read (a a p{x)) = — (M a p(x)), 

Now we replace the auxiliary fields cr uu (x) and cr dd (x) by spatially uniform condensates of the scalar mode as 
Cuu(^) = fui c<id(x) = &d- Then we find that cr u and ad correspond to the chiral condensates of the up and down 
quarks: a u = — (u a (x)u a (x)) /A c and o& = — (d a (x)d a (x)) /N c . We also replace a u d(x) and <Jdu(x) by spatially 
uniform condensates of the pseudoscalar mode as cr u d(a;) = ie(a;)7r, adu(x) — ie(x)Tr* . Then ir corresponds to the 
charged pion condensates: n = (iu a (cc)e(a;)d a (a;)) /A c and n* — (id a (x)e(x)u a (x)) /N c . 

Substituting these mean field values and Eqs. JjJJ, (|10|l in the action (|25p. we obtain 



S'[(p,X,X;^n,^d,^} = ^2 ~t~ V* +cr d + 2 kl 2 ) ^EEE X ( m ' f ) G 1 {m,n;x,(j) a )X{n,x) 



x a m,n 



where 



(26) 



X 



(m, a?) = (xL(m, x), d(m, x^J 



and 



X(n, x) 



u(n, x) 
Mn,x) 



(27) 



and 



G 1 (m,n;x,4> a ) 



M„ + isin k 



4>a(x) 



\e(x)-Tz8 ntm _ Nr / 2 



ie(x)-n*8 rhm _ NT / 2 



M d + i sin k 



4>a{x) 
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(28) 



M u = to + (d/2)er u and M d — to + {d/2)a d are the dynamical quark masses of up and down quarks. After performing 
the integration over the Grassmann variables X and X, we obtain the following effective free energy: 



F eff K,c7 d ,7r;T,/iB,Mi] = ~j~ K? + °d + 2 \A) ~T\ogl / dU TjDet^ 1 ^ 



Ay 



(29) 



where 



Det [G-\r 



)]-n 

n=l 



(l^j Kl 2 + { M u + isin (fc„ - iMi/ 2 )} ■ {Md - isin (k n + W 2 )} 



(30) 



with k n = k n + <f> a {x)/N T — i/iB/A c : Here we have rewritten the chemical potentials of the up and down quarks as 
Hb and fii using the definition /i u = hb/N c + Hi/2 and = Mb/Ac — (mj/2. The integration over Uq in the Polyakov 
gauge is defined in Appendix iBl 

It is impossible to perform the summation over n in Eq. I|30() analytically for general /iB and fii. However, we can 
do it in two interesting cases. One is the case that n\ is lower than a critical value fi\" ~ to^, where the pion does not 
condense, tt = fl5| . This case is relevant to realistic systems such as heavy ion collisions and electroneutral neutron 
stars. The other is the case that /ie = 0, which means we can put <t„ = cxd- In this case, we can compare our results 
with those obtained by the recent Monte Carlo lattice simulations Jlq. 
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B. The phase structure for /ii < 

First, we derive the analytical expression for the effective free energy in the case of [i\ < /i" 1 , where we can put 
7r = 0. The proof for the existence of such a critical value and its correspondence to the pion mass in the framework 
of strong coupling lattice QCD will be given in the next subsection. In this case, we can calculate the product over n 
using the formula in Appendix 1X1 



Det [G-h 



m, n; x, 



6 «)] = II { Mu + isin [ fc » + M^)/N T - i/i u ]} ■ {M d + isin [k n + (p a (x)/N T - i/i d ]} 



n=l 



(31) 



{2 cosh [N T E n ] + 2 cos[(j) a (x) - iN T (i u ]} • {2 cosh [N T E d ] + 2 cos[(j) a (x) - iN T fi d }} 



with one-dimensional up and down quark excitation energy E u — arcsinhAf u , E d — arcsinhAfj. 

Using the formula in Appendix 151 we can complete the SU(A r c ) integration over J7 in Eq. I|29[) for general N c 
Then we obtain the analytical expression for the effective free energy as follows: 



r rri . crii N c d , 2 

^eff[CT u ,(Td; J,/^B,MI < Ml J = (,°u 



Tlogl ^3 det Q 



n-\-i—j 



(32) 



where 

Qo = 4 cosh [E u /T] • cosh [E d /T] + 2 cosh [( Mu - f i d )/T] , 
Q ±1 = 2 cosh [EJT] • (cosh [/x d /T] ± sinh [fi d /T]) + 2 cosh [E d /T] • (cosh [^t u /T] ± sinh [fx u /T]) , (33) 
Q± 2 = cosh [(/i u + Hd)/T] ± sinh [(// u + ^ d )/T] , Q|r»|>3 = > 

and the determinant is to be taken with respect to i, j = 1, 2, . . . , iV c . 

Although the effective free energy has a complicated expression, it reduces to a simple form at T = 0. In the zero 
temperature limit N T — > oo, the contribution of 4> a {x) in Eq. I|31|) vanishes. Then the up quarks and the down quarks 
are completely decoupled, and we can write the effective free energy for T = as 



F eff [(7;T = 0,//B,Mi <Mi ri ] 



^pcr u - max {iV c /x u , A^ c £; u } 



a d - max { N c fi d , N c E d } 



(34) 



Each term coincides with the effective free energy at T = for JVf = 4 [Eq. (|18(l ] with quark chemical potential /i u or 
/id- Therefore the discussions on chiral restoration at finite density given in Sec. Ill Bl hold true in each sector. 



The phase diagrams 

At finite temperature, the integration over the temporal gauge link variable <p a of Eq. (|31|l nontrivially couples the 
up and down quark sectors. Let us see its consequences numerically here. We show in Fig. [3] the phase diagram in 
the plane of T and /j,b with N{ = 8, determined by minimizing the effective free energy Eq. (|32|) for N c — 3, d = 3 
with respect to the order parameters cr u , a d . 

First we consider the case of /ii = in which we can put er u = a d ; the phase diagram is shown in the left panel of 
Fig. |3 This should be compared to the case of four-flavor QCD in Fig. [21 The outline of the phase diagrams is quite 
consistent, and the discussions given in Sec. Ill Cl hold true here, while the critical temperature for Nf — 8 flavor is 
lower than that for Nf — 4. This is expected from a physical point of view: With more flavors, the chiral condensate 
is broken faster by their thermal excitations. In other words, the critical temperature decreases on increasing the 
number of quark flavors. What we want to emphasize here is that the integration over (j> a causes the N{ dependence 
of the critical temperature even at the mean field level. 

Next we consider the phase diagram at [i\ = 0.2 with m — 0.4. As shown in the right panel of Fig- El the introduction 
of fii splits the first order phase transition line into two lines associated with the jump of the chiral condensate of 
the up quark and that of the down quark. They shift in opposite directions of /j,b and terminate at the second 
order critical end points. Such a successive chiral restoration with finite /ii on increasing /ie is easy to understand 
intuitively from the fact that up and down quarks experience the chemical potential differently, /i u = /j,b/N c + /Xi/2 
and /x d = hb/N c - 

What is nontrivial in our result is the position of the critical end points. The temperatures of the critical end points 
for er u and a d are both higher than that in the case of fii = 0. This is in contrast to studies using the random matrix 
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FIG. 3: Phase diagrams of strong coupling lattice QCD with N c = 3, N{ = 8 in the T-/^b plane with zero (the left panel) and 
finite (the right one) isospin chemical potential fix. In the left panel, the solid (dashed) line denotes the phase transition line 
for the first (second) order phase transition in the chiral limit at fj,x = 0. The first order line and its critical end point for 
m = 0.4 are also shown. The right panel shows the phase diagram at fix = 0.2 and m = 0.4. The two solid lines represent 
the phase transition lines for the first order phase transitions associated with up (the left line) and down (the right one) quark 
condensates, respectively. The first order phase transition line for fix = is plotted by the dotted line for comparison. The 
symbol <r u or <7d in the figure represents the region where the condensate has a large value. 



model [35j] , Nambu— Jona-Lasinio model j3|| and ladder QCD approach (3?|], where the temperature of the critical 
end points is not affected by fi\ because the up quark sector and the down quark sector are completely uncoupled 
in such approaches. In our approach, as discussed above, the up and down quarks are nontrivially coupled via the 
integration of the temporal gauge link variable Uq at finite temperature. 



C. The phase structure for fin — 



We derive the analytical expression for the effective free energy in the case of /ib = 0. In this case, we can put 
a = <t u = <7d and obtain the following expression for Eq. I|30|l : 



Det [G^(m,n;x,cb a )] = \ \ 



d Y I |2 , far , ■ • (1 , f nr • • (, , <t>a{x) . .Ml 

— 7r + < ill + 1 sin k n H 1 — > • < M — 1 sin k n -\ + 1 — 



(35) 

This expression has the same form as the determinant calculated in [2^ for SU C (2) QCD. Therefore we can use the 
formula in the Appendix of 23] for the product over the Matsubara frequencies, and then we obtain 

Det [G^im^nix^a)] = {2 cosh [N T E-] + 2 cos [(/> a (x)]} ■ {2 cosh [N T E+] + 2 cos [<p a {x)]} (36) 

with one-dimensional up quark and down antiquark (or up antiquark and down quark) excitation energy 



E± =arccosh^(l + Af 2 )cosh^[ W /2] + (d/2) 2 |7r| 2 ±Msinh[ W /2]J . (37) 

Using the formula in Appendix iBl we can complete the SU(iV c ) integration over U$ in Eq. I|29|l for general N c . 
Then we obtain the analytical expression for the effective free energy as follows: 



F eS [a, tt; T, = 0, M i] - ^ (cr 2 + |tt| 2 ) - Tlog |^ det R n+l - 3 1 



(38) 



where 

R Q = 4 cosh [E-/T] ■ cosh [E+/T] + 2 , 
R±! = 2cosh[£;_/r] + 2cosh[£; + /T] , (39) 

R±2 = 1 , R\n\>3 = , 
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and the determinant is to be taken with respect to i, j = 1, 2, . . . , N c . Note that for m — and p\ = 0, the effective 
free energy is a function of only a 2 + \ir\ 2 . Of course, this is from the chiral symmetry U(2)a of the original action at 
m = p,\ = 0. The symmetry between a and it will play an important role in understanding numerical results on the 
phase structure. 

Although the effective free energy has a complicated expression, it reduces to a simple form at T = 0. In the zero 
temperature limit N T —>■ oo, the contribution of 4> a {x) in Eq. 1(36(1 vanishes, and we obtain the effective free energy 
for T = as 



F cS [a, 7r; T = 0, p B = 0, pi] = N c 



- (E- + E+) 



(40) 



If we replace it by A and p\ by pb, the above expression corresponds to the free energy at T — of SU C (2) QCD 
studied in |23j ] up to the overall factor N c . This correspondence is rather natural from the point of view of the 
symmetry and the effect of the chemical potential. In the present system, up and down quarks arc indistinguishable 
because of chiral symmetry at p\ = and they undergo the effect of p\ with opposite signs. On the other hand, in 
SU C (2) QCD, the quark and antiquark are indistinguishable because of the Pauli-Gursey symmetry ^S$^ at pb = 0, 
and they also undergo the effect of pb with opposite signs. Therefore there is a correspondence between the up (down) 
quark in the finite isospin density system and the (anti)quark in the finite baryon density SU C (2) system, which also 
means a correspondence between the pion condensation and the diquark condensation. Note that the difference of the 
effective free energy in the two systems at T ^ is due to the integration of Uq for the finite isospin density system, 
which couples up and down quarks nontrivially at finite temperature. 

As calculated in [23( for the case of SU C (2) QCD, we can derive analytical properties from the effective free energy 
at T = in Eq. ((40(1 . The system has two critical chemical potentials, the lower one p l ° w and the upper one p^ p , 
given, respectively, by 



low 



= 2 arccoshV 1 + mM and /i" p = 2 arccoshV 1 + K 2 , 



(41) 



where M is a solution of the chiral gap equation with quark mass m, and we have defined K as the solution of the 
equation 



K 



m 
~K 



a + K 2 



-1/2 



(42) 



low 



The nonvanishing value of the pion condensate 7r ^ is possible only for p x ° w < /jj < /i" p . As long as < fi' c 
the empty vacuum gives p\ = —dF c s /dfii = and n = 0. On the other hand, for fii > /i" p , saturation of the isospin 
density occurs, leading to p\ — N c and n = 0. For sufficiently small m, the lower critical chemical potential reduces 
to 



low 



= 2m 1 / 2 



{l + d 2 ) 1 ' 2 -! 



1/4 



(43) 



This expression can be rewritten as 



low 



ran with being the pion mass obtained from the excitation spectrum 



in the vacuum 28] up to leading order of the 1/d expansion. This is nothing but a proof of the fact that the critical 
value of the isospin chemical potential n" 1 assumed in the previous subsection indeed exists and corresponds to the 
pion mass. 



Condensates and isospin density 



Now we determine the chiral condensate a and the pion condensate 7r numerically by minimizing the effective free 
energy with Nf = 8 flavors in Eq. (|38|l for iV c — 3, d — 3. The isospin density p\ = —dF c s/dp\ is also calculated as a 
function of pi and T. The results are shown in Fig. 0] for m = and m — 0.02. Here we choose it to be real. 

Let us first consider a, ir and p\ as functions of p.\ at T = 0, shown in the left panels of Fig. In the chiral limit 
772 = (the upper left panel), a is always zero, while 7r decreases monotonically as a function of pi and shows a second 
order transition when p\ becomes of order 2. On the other hand, p\ increases until the saturation point p\ = 3 where 
quarks occupy the maximally allowed configurations by the Fermi statistics. 

The lower left panel of Fig. 0] is the result for small quark mass m = 0.02. As we discussed above, there exists 
a lower critical chemical potential /jJ. ow ~ given by Eq. 1)41(1. Both 7r and p\ start to take finite values only 
for m > p}° w . One can view the behavior of a and n with finite quark mass as the manifestation of two different 
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FIG. 4: Chiral condensate a, pion condensate n and isospin density pi for N c = 3, Nt = 8. In the left panels, they are plotted 
as functions of the isospin chemical potential /xi at T = for m = (the upper panel) and m = 0.02 (the lower one). In the 
right panels, they are plotted as functions of the temperature T at zxi — 0.4 for m = (the upper panel) and m = 0.02 (the 
lower one). The dotted lines in the lower panels indicate the total magnitude of the condensates \/ a 2 + ty 2 . 



mechanisms. One is a continuous "rotation" from chiral condensation to pion condensation above /xi = /x c ow with 
\/ o 2 + 7r 2 varying smoothly. The other is the "saturation effect" ; the isospin density forces the pion condensate to 
decrease and disappear for large /xt as seen in the previous case of m = 0. 

The "rotation" can be understood as follows. As we have discussed, the effective free energy at vanishing m and 
/ii has a chiral symmetry between the chiral condensate and the pion condensate. The effect of m (/ii) is to break 
this symmetry in the direction of the chiral (pion) condensation favored. Therefore at finite m a relatively large a 
appears predominantly for the small /ii region. Once fXi exceeds a threshold value /xj° w , a decreases while 7r increases, 
because the effect of /ii surpasses that of m. As /ii becomes larger, the pion condensate begins to decrease in turn by 
the effect of the "saturation" and eventually disappears when /xi exceeds the upper critical value /x" p (of the order of 
2 for T — 0) . Such a behavior of 7r is also observed in recent Monte Carlo simulations of the same system . 

Next we consider the chiral and pion condensates as functions of T at /ii = 0.4 in the right panels of Fig. In 
the chiral limit m = (the upper right panel), a is always zero, while 7r decreases monotonically and shows a second 
order phase transition at some critical temperature T c . 

For small quark mass m = 0.02 (the lower right panel), both a and it have finite values. 7r decreases monotonically 
as T increases and shows a second order phase transition at a critical temperature smaller than for m = 0. On the 
other hand, a increases as 7r decreases, so that the total condensate \/ a 2 + it 2 is a smoothly varying function of T. 
The understanding based on the approximate symmetry between the chiral and pion condensates is thus valid. An 
interesting observation is that a, although it is a continuous function of T, has a cusp shape associated with the phase 
transition of the pion condensate. 

Finally let us compare the results for m = and those for m = 0.02 in Fig. Looking at the two figures at the 
same T or /xi, we find that the pion condensate it for m — and the total condensate \/ o 2 + tx 2 for rn = 0.02 have 
almost the same behavior. This indicates that, although the current quark mass suppresses 7r, the price to pay is an 
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FIG. 5: Phase diagram of strong coupling lattice QCD with 
iV c = 3, Nf = 8 in the T-/j,i plane. The solid (dashed) line 
denotes the second order critical line for pion condensation 
for m = (m = 0.02). 
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FIG. 6: Phase structure of strong coupling lattice QCD with 
N c = 3, Nf = 8 in the T-[i\-m space. The surface represents 
the critical surface for pion condensation, which separates 
the region where n =^ (inside) from the region where 7r = 
(outside) . 



increase in a, so as to make the total condensate insensitive to the presence of small to. 



The phase diagrams 

Now we show in Fig.Elthe phase diagram of the strong coupling lattice QCD with N c = 3, Nf = 8 in the T-fi\ plane. 
The solid line denotes the critical line separating the pion condensation phase it ^ from the normal phase it = in 
the chiral limit to — and the dashed line is for small quark mass m = 0.02. The phase transition is of second order 
on these critical lines. The chiral condensate a is everywhere zero for m = 0, while it is everywhere finite for m^0. 
In the latter case, however, a is particularly large in the region between the solid line and the dashed line. 

Then the phase structure in the three dimensional T-jii-m space is shown in Fig. The pion condensate has a 
nonvanishing value inside the critical surface and the phase transition is of second order everywhere on this critical 
surface. The second order phase transition is consistent with other analyses employing the mean field approximation 
[35ll3^| . On the other hand, Monte Carlo simulations of QCD with finite fii [ly show an indication of the presence 
of a tricritical point in the T-jii plane at which the property of the phase transition line changes from second order to 
first order as p\ increases. Clarifying what induces such a first order phase transition is an interesting open problem. 

The plane at T — in Fig.|S|shows the phase diagram in terms of p\ and to. The lower left of the plane corresponds 
to the vacuum with no isospin density present, p\ = 0. On the other hand, the upper right of the plane corresponds 
to the saturated system, p\ = 3, in which every lattice site is occupied by three pions. There is a pion condensation 
phase (a superfluid phase of isospin charge) with ir ^ and < p\ < 3 bounded by the above two limiting cases. 
They are separated by the two critical lines given in Eq. 1)41(1. 

It is worth mentioning here that we can see the corresponding system in the context of condensed matter physics; 
the hard-core boson Hubbard model has a similar phase diagram in which a superfluid phase is sandwiched by Mott- 
insulating phases with zero or full density |40| . Here we can give a physical explanation for the saturation effect, 
that is, why the saturated isospin density forces the pion condensate to disappear. Suppose that one imposes a small 
external field for the isospin charge on the system in which every lattice site is occupied maximally by pions. However, 
a pion at one site cannot hop to the next site, owing to the Pauli principle of constituent quarks, and thus the isospin 
current never appears. This results in a zero condensate because the superfluid current is proportional to the square 
of the absolute value of the condensate. This is nothing but an analogous phenomenon to the Mott insulator. 

The fact that SU C (3) QCD at finite isospin density has the same phase structure as that of SU C (2) QCD at finite 
baryon density can be understood from the aspect of the symmetry and its breaking pattern. SU C (3) QCD at pi = 
and to = has chiral symmetry between the chiral condensate and the pion condensate. On the other hand, SU C (2) 
QCD at /ib = and to = has Pauli-Gursey symmetry between the chiral condensate and the diquark condensate. 
In both theories, the introduction of m breaks the symmetries so that the chiral condensate is favored, and the 
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introduction of /xi//zb breaks the symmetries so that the pion/diquark condensate is favored. This is the reason why 
these two theories have similar phase structures at finite density. 

IV. SUMMARY AND DISCUSSION 

In this paper, we studied the phase structure of hot and dense QCD with N c — 3 and Nf = 4, 8 in the strong 
coupling limit, formulated on a lattice. We employed the 1/d expansion (only in the spatial direction) and the mean 
field approximation to derive the free energy written in terms of the chiral condensate a and the pion condensate it 
at finite temperature T, baryon chemical potential /iB, isospin chemical potential fii and quark mass m. 

In Sec. [n] we investigated the phase structure with N{ — 4 in the T-fiB plane, and the phase transition of chiral 
restoration was found to be of second order in the high temperature region and of first order in the low temperature 
region. Analytical formulas for the critical line of the second order transition and the position of the tricritical 
point were derived. The critical temperature at small quark chemical potential fi — hb/N c can be expanded as 
T c (n) ~ T c (0) - lW/TJO), and the slope was found to be much larger than that calculated in Monte Carlo lattice 
simulations 0, El El E| ■ We discussed how this difference could be understood as a result of the strong coupling 
limit (j — ► oo and the chiral limit m = 0. We also studied the effect of finite quark mass on the phase diagram: 
The position of the critical end point has a relatively large dependence on the quark mass. In paticular for small to, 
the critical end point shifts in the direction of smaller T and smaller /iB. This is because the gradient of the phase 
transition line near the tricritical point is positive; how this feature in the phase diagram of QCD is affected by the 
introduction of dynamical baryons beyond the leading order of the 1/d expansion should be studied in future work. 

In Sec. IIIII we derived the free energy including /ii by extending the formulation studied in Sec. ITU with two species 
of staggered fermion (Nf — 8). First, we observed that the critical temperature of chiral restoration decreases on 
increasing the number of quark flavors due to the thermal excitations of mesons. Then we studied the effect of /zi on 
the phase diagram in the T-fiB plane at finite quark mass: The introduction of fi\ splits and shifts the first order phase 
transition line in opposite directions of /iB- In addition, /ij moves the two critical end points to higher temperatures 
than that for /ii = 0. This is a nontrivial effect originating from integration over the temporal gauge link variable at 
finite temperature, which couples the up and down quark sectors. 

Finally, we studied the phase structure in the space of T, fii and m and found a formal correspondence between color 
SU C (3) QCD with finite isospin density and color SU C (2) QCD with finite baryon density. Such a correspondence can 
be understood by considering the symmetry and its breaking pattern: chiral (Pauli-Gursey) symmetry at to = and 
/ii = (^b = 0) is broken by the quark mass in the direction of chiral condensation, while it is broken by the isospin 
(baryon) chemical potential in the direction of pion (diquark) condensation. Although our results are in principle 
limited to strong coupling, the behavior of a, w, and p\ in the T-^\-m space has qualitative agreement with recent 
lattice data. We also found that the disappearance of pion condensation with saturated density can be understood 
as a Mott-insulating phenomenon. 

The remaining work to be explored is to investigate the phase structure with finite /ib and finite fii as well as with T 
and m, where a nontrivial competition between the chiral and pion condensates takes place. In our formulation, such 
a study is possible by using the effective free energy given in Eq. (|29|l with numerical summation over the Matsubara 
frequency and numerical integration over the temporal gauge link variable. How the chiral condensate at finite (Ib is 
replaced by the pion condensate as jii increases is one of the most interesting problem to be examined. 
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APPENDIX A: SUMMATION OVER THE MATSUBARA FREQUENCY 

In this appendix, we give the formula for the product over the Matsubara frequencies of the expression 

N T /2 N T 

I = [sin 2 k n + A 2 ] = J] [sin 2 k n + A 2 ] 1/2 (Al) 

n— 1 n— 1 
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with k n = 2ir(n — 1/2) /N T + (f>/N T — ifj, and N T an even integer. Let us first take the logarithm of the above expression 
and differentiate it with respect to A: 



d_ 



log/=^- 



A 



n=l 



sin fc„ + A 2 



(A2) 



Because sin(fc„ + 4>/N T — i/x) is invariant under the shift n — > rt + N T , we can make the summation of n over the range 
n — — oo to n = +00 with an appropriate degeneracy factor f2. Then the residue theorem enables us to replace the 
summation by a complex integral as 



1 T 1 



A -iN T ^ 

k I at ;..\ 1 \2\ + c iN T z 2-~i 



A 



2ni sin 2 (z + 4>/N T - i/x) + A 2 



2 sin(z + 0/A T - ifj.) cos(z + 0/7V T - i/x) 1 + e iAr - ? 



(A3) 



Owing to the infinite range of the summation over n, we can choose the closed contour at infinity for the complex 
integrals with respect to z, and thus such complex integrals go to zero, z are the residues satisfying 

sin 2 (z + (j>/N T - ifi) + A 2 = . 



Solving this equation, we obtain 

z + <j>/N T - in = ±\E + 2-kti , ±iE + it + 2im , 
with E — arcsinhA and n = — 00, . . . , 00. Substituting Eq. I|A5(I into Eq. (|A3|) . we obtain 



(A4) 
(A5) 



^A ° g O rnH V 



£1 ^— ' cosh E 

n— — 00 



as 

9A 
_9_ 

dA 

_a_ 

dA 



I _|_ e -N T E-i<l>-N T ^ 1 -f e N T E-i<p-N Tt i 



I _|_ Q — N T E — i{p—N T [j. 1 _|_ e N T E-i<f)-N Tt i 



(A6) 



[log{ 



AS _L „-i0-Af T M 



} + log {< 



-JV T B , e -i0-JV T/ i 



}] 



(-i0 - A^ T ^ + log {2 cosh LY T .E] + 2 cos[0 - iiV Tj u]}) 



We note that the degeneracy factor Q is just canceled by the infinite degeneracy of the summation on n. After 
integration with respect to A, we find that the result is expressed in a rather simple form up to irrelevant factors, 



1 = 2 cosh [N T E] + 2 cos[0 - iN T fx] , 



(A7) 



with E = arcsinhA. 



APPENDIX B: SU(N C ) INTEGRATION IN THE POLYAKOV GAUGE 

In this appendix, we give the formula for the SU(iV c ) integration 

J = JdU F(U ) . 

Here we consider the function that can be written as the product of f(6 a ) in the Polyakov gauge as 

F(U o ) = l[f(0 a ) . 

a=l 

Since the SU(iV c ) Haar measure in the Polyakov gauge is given by 

J J -*a=l Z7T a<b 



7T J 1 

2tt 



-,i<j>i{N c -ii) . . . p it/>N c (N c -i N€ ) 



■ 2 



(Bl) 



(B2) 



(B3) 
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the integration becomes 



J 



n 



2tt 



a i0i(iVc-u) . . . p i0N c (Wc-ijv c ) 



2^ x[]/(« a ) 

\a=l / a=l 



Because the delta function can be rewritten as 



(N c \ oo 
^2 <t>a \ = exp 
a—1 / n— — oo 



-in ^ (?!> a 

a=l 



we obtain 



uo i\ c ^ UO 

n— — oo a—1 n— — oo 

Here we have defined M„ as 

M„= / ^/(0)exp[-in0] , 
and the determinant is to be taken with respect to i,j = 1, 2, . . . , iV c . 



(B4) 



(B5) 



(B6) 



(B7) 
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